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In this proceeding, we report our recent study of the mixing mechanism between D s \ (2460) and 
D s \ (2536). On the basis of Godfrey-Isgur model, we consider the D*K hadron loop effect on the 
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we can extract the poles as well as the mixing angles. Through this method, we simultaneously 
determine the masses, widths and mixing angle of these two physical states and the results agree 
well with the experimental measurement. Besides the mass shift, we also find that the hadron 
loop effects can cause a significant shift for the mixing angles from the value determined in the 
heavy quark symmetry limit. 
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1. Introduction 

Last year, BaBar Collaboration updated their measurement of the mass and width of D s \ (2536) 
meson [1]. Before this measurement, only an up limit of the width had been set [2]. This new 
measurement should be useful for the study of the internal structures of D 5 i(2536) and its axial- 
vector partner D s \ (2460). Thus, it initiates a lot of theoretical interests in the mixing of these two 
states. 

We recall first some observations concerning the P-wave D s state. Theoretically, the Godfrey - 
Isgur (GI) model [3] has made a great success in providing an overall description of charmonium 
and charmed meson spectra. However, it is also realized that there exist apparent discrepancies 
between theory and experiment when a state is located close to an open threshold. In the D s 
spectrum, the lowest 5-wave states D s and D* can be well explained by the GI model. But for those 
four P-wave states, i.e. D, (2317), Ai(2460), D,i(2536), and D l2 (2573), there are significant 
discrepancies between the quark model prediction and experimental observation. For instance, the 
masses of D 4 -o(2317) and D s i(2460) are about 100 MeV lower than the quark model prediction, 
and a dynamical reason for such large mass shifts may be due to the DK and D*K thresholds, 
respectively [4, 5, 6, 7]. 

State mixing makes these two axial- vector states D s \ (2460) and D s \ (2536) more interesting. 
Such a heavy-light quark-antiquark system could be an ideal place to test the heavy quark symme- 
try. Namely, in the heavy quark symmetry limit, D v i(2460) is a pure j = 1/2 state and couples 
to D*K channel through an S wave, while D s \ (2536) is a pure 7 = 3/2 state and couples to D*K 
channel through a D wave. The heavy quark symmetry is expected to be broken at the order of 
l/m c . Explicit introduction of a broken mechanism for the heavy quark symmetry in these two 
axial-vector states is thus necessary. 

Since both D s \ (2460) and D s \ (2536) are not charge conjugation eigenstates, and they both 
have strong 5-wave couplings \.oD*K, D* s r\, and DK*, we can then assume D 5 \ (2460) andD 5 i (2536) 
to be mixed states of the pure 3 P\ and l P\ cs states through intermediate meson loops. We mention 
that the D*K channel plays a dominant role in such a mechanism. By calculating the intermediate 
meson loop transitions, we can simultaneously determine the masses, widths and mixing angles of 
the two physical states. A similar mechanism has been applied to the ao(980) — /o(980) mixing [8]. 

We mention that many other scenarios are also proposed to explain the experimental data for 
D 4 i(2460) and D s i(2536), such as D*K molecule, tetra-quark state and dynamically generated 
state [9, 10, 11, 12]. Besides the mass shift, the hadron loop effect as a general mechanism can 
also affect the decay and production process. For instance, it has been used to evade the helicity 
selection rule in % c \ — > VV and Xc2 — > VP [13], account for the large branch ratio of the non-DD 
decay of i//"(3770) [14, 15], explain the large isospin violation process T](1405/1475) — > 3n [16], 
and predict a direct measurement of open threshold effects in e + e~ — > J/ [17]. 

2. Mixing through intermediate meson loops 

In this part, we show our basic formulas for the mixing scheme through intermediate meson 
loops. If two states \d) and \b) can transit to each other through hadron loops, the propagator of 
such a two-state system can be written as a 2 x 2 matrix G a t [18]. In general, the propagator matrix 
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has two poles in the complex energy plane, which correspond to the physical states. From the pole 
position, we can determine the masses and widths of the physical states. The physical states \A) 
and \B) should be a mixture of \a) and \b), 

=R(e,<t>)(Yl\ (2.i) 





\b) / 

where 6 is the mixing angle, and <j> is a possible relative phase between \a) and \b). Then the 
propagator matrix of \A) and \B) becomes Gab = RG ab PJ and should be diagonal. So we can 
obtain the mixing parameters (9, </>) by diagonalizing G ab . 

Here, we show the propagator matrix. For the spin-0 states, i.e. scalar or pseudo-scalar, the 
propagator matrix is 

G ab = l —( DbDab ) , (2.2) 

D a D b -D 2 ab \D ba D a ) 

where D a and D b are the denominators of the single propagators of \a) and \b), respectively, and 
D ab is the mixing term. For the spin-1 states, i.e. vector and axial- vector, the propagator matrix is 

G ^~ lP d^0M +lQ detG^' (23) 

where P^ v = g^ v — p^p v jp 2 and Q^ v = p^p v / p 2 are the transverse and longitudinal projector, 
respectively. The poles and mixing angles are only related to the transverse part and 



,2 _ x ( m 2 b + Yl b (s)-s -U ab (s) 

-U ab (s) m 2 a + U a (s)-s 



G ab (s) =M£ b -5 ab s= | ^ 2 , tt r.\ _ I ' (2.4) 



where Yl a , Yl b and Yl ab are the transverse coefficients of the self-energy functions and the mixing 
term, respectively. 

For the case of D s \ (2460) and D s \ (2536), we apply the following mixing scheme 
10,1(2460)) = cosepPO-sinee^Pi) 

|fl,i(2536)) = smde-^^P^+cose^Pi) . (2.5) 

In the heavy quark limit, the ideal mixing angle is 6q = 35.26° in our convention. Considering par- 
ity conservation and the OZI rule, the important intermediate states that can couple to D s \ (2460) 
and D s i(2536) are D*K, D*r\ and DK* . To compute these diagrams, we need to know the cou- 
plings of the vertices as well as form factors to remove the ultraviolet (UV) divergences in the loop 
integrals. 

3. Result and discussion 

At hadronic level, the Axial-vector- Vector-Pseudoscalar vertex can be written as i(gsg >lv + 
goP fl p v ), where g$ and go are the S and Z)-wave couplings. Near threshold, the D-wave couplings 
are suppressed due to the small momentum. So we only consider the 5-wave couplings, and use 
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the chiral quark model to evaluate g$ [19, 20, 21, 22]. In the chiral quark model, the light pseudo- 
scalars and light vectors are treated as chiral fields. 

Our result for the S- wave couplings are shown in Fig. 1 . We can see that they are insensitive to 
the initial meson mass. The couplings to D*K are strong, the couplings to D*r\ are rather weak, and 
the couplings to DK* are vanishing as a leading approximation. To remove the UV divergences in 
the loop integrals, we modify the non-relativistic (NR.) exponential form factors in the quark model 
to a covariant form, i.e. exp(— q 2 /4a 2 ) — > exp(g 2 — m 2 /A 2 ), where A is the cutoff energy and can 
be fixed by the chiral quark model. 




Figure 1: The absolute values of coupling g s as Figure 2: The mixing term U ab with different 
a function of the initial meson mass M,- [18]. cutoff values for the exponential form factor. 

In the propagator matrix Eqs. (2.3-2.4), what we need to know are the self-energy functions IT, 
and Tib, and the mixing term Yl a b- They are found to be proportional to each other [18]. In Fig. 2 
we show the energy dependence of Yi a b- The blue lines are the results with cutoff A = 1.174 GeV 
fixed by the chiral quark model. We can see that both the real and the imaginary part have two kinks 
due to the charged and neutral D*K thresholds. Below the thresholds the imaginary part is zero, 
while above the thresholds it increases quickly. At the energy near 2460 MeV, only the real part 
contributes. At the energy near 2536 MeV, both the real and imaginary parts have contributions. 
We also vary A to see the cutoff dependence and find that the real part is sensitive to A, while the 
imaginary part is unchanged. 
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Figure 3: Pole structures highlighted by the zero Figure 4: Schematic plot for the mass-shift pro- 
values of det[G] in the propagator matrix. cedure [18]. 

The poles in the propagator matrix G^ b v are equal to the zero points in G a b as shown in Fig. 3. 
For the bare cs masses, we adopt the values from the GI model, m[ 3 Pi] = 2.57 GeV and mpPi] = 
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2.53 GeV [3]. From Fig. 3 we can see two poles: one is near 2460 MeV and the other is near 
2536 MeV. The higher pole is insensitive to the cutoff while the lower pole turns out to be rather 
sensitive. 

We list the masses and widths in Table 1 with A = 1.174 GeV. The results agree with ex- 
periment, which indicate that the exponential form factor from the quark model can give a good 
estimate of the real part. The detailed mass shift procedures of these two physical states are demon- 
strated in Fig. 4. 

Table 1: Masses and widths obtained from the pole analysis [18]. 



[m - /£] (MeV) 


Ai(2460) 


Avi(2536) 


Kb 


2454.5 


2544.9-1.0/ 


K b +K b 


2455.8 


2544.9-1.1/ 


Experiment 


2459.5 


2535.08-0.46/ 



Besides the masses and widths, we can also extract the mixing parameters and the results are 
shown in Table 2. We can see that the mixing angles are different for these two physical states at 
their own mass shells. It is also interesting to see the deviation of the mixing angles from the heavy 
quark limit. For D v i(2460) the deviation is large A9 = 12.3°, and for D s i(2536) the deviation is 
relatively small A9 = 4.4°. The new experimental data from BaBar put a strong constraint on the 
mixing angle of D v i(2536), but with two solutions, i.e. B\ = 32.1° or 62 = 38.4° [18], which are 
symmetric to the ideal mixing angle. It shows that our theoretical analysis favors the bigger one. 



Table 2: The mixing angle 8 and relative phase (j) extracted at the two poles in those three diagonalization 
schemes [18]. 





Ai(2460) 


Ai(2536) 


I 


II 


III 


I 


II 


III 




{47.5, 0} 


{47.5, 0} 


{47.5, 0} 


{39.7, -6.4} 


{39.7, 6.4} 


{39.7, 0} 


Kb+K, 


{47.6, 0} 


{47.6, 0} 


{47.6, 0} 


{39.8, -6.5} 


{39.8, 6.5} 


{39.7, 0} 



4. Summary 

When taking into account the D*K loop corrections in the Godfrey-Isgur model, we can ex- 
plain the masses and widths of the physical states D s i(2460) and Z) s [(2536), and extract their 
mixing angles. We also find that when there are strong 5-wave coupling channels, the hadron loop 
corrections can cause both large mass shifts from the quark model prediction and large mixing 
angle shifts from the heavy quark limit. Also, from our calculation we can see that the exponential 
form factor from the quark model can give a good estimate of the real part of the meson loops. 
Our study may provide further insights into the role of hadron loops as an important unquenched 
correction to the potential quark model. 
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